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O Abstract 
(N 

^ We have performed a theoretical study of electronic transport in single and bilayer graphene based on the standard linear-response 
(Kubo) formalism and continuum-model descriptions of the graphene band structure. We are focusing especially on the interband 
^"5 contribution to the optical conductivity cr(w). Analytical results are obtained for a variety of situations, which allow clear identifi- 
0^ pation of features in cr(w) that are associated with relevant electronic energy scales. Our work extends previous numerical studies 
and elucidates ways to infer electronic properties of graphene samples from optical-conductivity measurements. 
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1. Introduction 

. Graphene is a single sheet of carbon atoms tightly packed 
Into a two-dimensional (2D) honeycomb lattice. This mate- 
rial has recently become available for experimental study 1 1, 2], 
and its exotic electronic properties are attracting a lot of inter- 
est f^. In particular, the conical shape of conduction and va- 
lence bands, together with the absence of a gap, near the K and 
K' points in the Brillouin zone renders graphene an intriguing 
type of quasi-relativistic condensed-matter system |4]. Recent 
experiments have verified that the band dispersion of charge 
carriers in graphene is indeed linear as expected for massless 
Dirac fermions jslSl- A multitude of interesting physical ef- 
fects arising in single-layer and bilayer graphene samples have 
been discussed |^, both theoretically and experimentally. 

In this work, we focus on the ac electric transport properties 
of graphene, which have been the subject of numerous theo- 
retical (mostly numerical) studies Il7l 428tl and several recent ex- 
periments Measurements of quantities related to the 
optical conductivity are expected to give deeper insight into the 
electronic properties of graphene samples, making it possible to 
infer details of their morphology [ 28l 13311 and suitability for ap- 
plications. It is thus important to obtain a clear understanding of 
the features exhibited in the frequency-dependent conductivity, 
in particular, their relation to microscopic parameters and be- 
havior at finite temperature T . Furthermore, it is advantageous 
to have mathematical expressions available that can be straight- 
forwardly used for comparison with data. This is the motivation 
for our study. We have developed a formalism that lends itself 
for generalization to many situations, in particular, the treat- 
ment of inelastic scattering. Extrapolation to zero frequency 
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will enable us to discuss the dc conductivity, shedding new light 
on the phenomenon of minimal conductivity in graphene. The 
full range of our results will be presented elsewhere; here we 
focus on discussing the method and present selected results. 

2. Calculation of the conductivity: Basic formalism 

Our starting point is the familiar Issll Kubo formula 



cr^v(w) 



yJ-C 



dt e 



/(m-;0+)/ 



K. 



with the kernel 



V= l-Tr{e-'-;^e'-[r,,p]} 



(1) 



(2) 



Here - -er^ = i[H , r^] is the current operator, and p the 
density matrix. An alternative expression for the kernel ll35ll . 



K^y = £''dA Tr{e-'f e'^" g 



(3) 



will become particularly useful to enable discussion of the ef- 
fect of inelastic scattering. When quasiparticle interactions 
are neglected, it is possible f]3\ to express the conductivity in 
terms of matrix elements between eigenstates | n> of the single- 
particle Hamiltonian having energy e„: 



V (n\[H,r,]\n'){n'\[H,ry]\n) 

n,n' 



X [/(£„) -/(£„')]. (4) 



Here /(e) = 1/(1-1- exp{[e - ^i\l{k^T]} denotes the Fermi func- 
tion, which depends on the chemical potential yu. In the follow- 
ing, we use the expression (|4]i to derive conductivity formulae 
applicable to graphene. 
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2.1. Clean limit: Plane-wave representation 

Using continuum-model descriptions of the band structure 
near the K-points, single-particle eigenstates of clean graphene 
systems can be written as a direct product of a plane wave in 
real space and a 2A^-spinor (the latter subsuming the two sub- 
lattice and layer-index degrees of freedom): \ n) - | k) (g) | <j\. 
In single-layer graphene, cr distinguishes the two (electron and 
hole) bands. Note that the spinor wave function depends on 
wave vector k. The current operator y'^ or, equivalently, the 
commutator [H , is diagonal in the real-space part | k) while 
having off-diagonal matrix elements in pseudospin space. It is 
then straightforward to specialize Eq. (|4]i to this case, finding 



2 2N „ ,2, 

e r a I 



cPk <cr|w^(k)|(r')k<tr'|Wv(k)|(r)k 



■)2 fiko-' - fiktr + h(jL)- /0+ 



^/(gktr) -/(gko-O 

S 



(5) 



Here H'^(k) is found from {H , r^] |k) - vv^(k)|k). Two con- 
tributions to the conductivity can be distinguished, arising from 
terms with cr - cr' and cr + cr' , respectively. It is customary to 
call these the ;nfra-band and /nfer-band contributions. Defining 
w;^-'(k) = (cr|w^(k)|cr')k, we find 



cr, 



(intra) 



O"0 



(inter) 



2ha) 



2/' 



^ (k)wrHk)/'(ek<r), 
(6) 

d^k S (hcj - [ek<T - fiktr']) 



-w<r"''(k) wr'(k) 



■cosh(5|^) + cosh(2-i^)' 



(7) 



For brevity, we use the scale factor ctq - ge^/(2nh), where 
g = 4 has been introduced to account for the quasiparticle de- 
generacy (real spin and valley) in graphene. The intra-band 
term (|6]l is the usual dc Drude conductivity, which depends on 
states in the vicinity of the Fermi surface where the derivative 
/' of the Fermi function is peaked. It vanishes at the neutral- 
ity point (fi = 0) in the zero-temperature limit. The inter-band 
contribution d?) is the interesting part for finite oj. It is cal- 
culated straightforwardly using continuum-model descriptions 
of single-layer and bilayer graphene. The expression (|7]i given 
here is well-suited for obtaining analytical results for the de- 
pendence on temperature but cannot be used to go beyond the 
clean limit. To discuss the effect of disorder, a more general 
formula is needed that will be given in the next subsection. 

2.2. General conductivity formula in terms of Greens functions 

Mathematical manipulation of Eq. (|4]l yields the conductivity 
expressed in terms of single-particle Greens functions 
I22I1 . It has the general form 



luctiyity 

Elil 



g-^y(^) 



de Tfjvie, w) 



f\s^--f(s-- 



(8) 



We have derived a new and, for our purposes, more convenient 
expression for the diagonal part of the frequency-dependent 
transmission function, 

(hcij , htij 

— -s-iO\-—-s-iO' 



-F,\—-s-iQ' 

+C.C. , 



■ E + iO' 



given here in terms of functions 



(9) 



(10) 



G(r, r';z) = G(r-r',0;z) -: Gr-r'(z) is the real-space represen- 
tation of the single-particle Greens function in a translationally 
invariant system, and the trace Trsi is performed only over sub- 
lattice and layer degrees of freedom. 

It is straightforward to specialize the general conductivity 
formula obtained in this subsection to the clean limit. Perform- 
ing a Fourier transformation and using the fact that the single- 
particle Hamiltonian // — > //k becomes diagonal in real space, 
we find 



Fi,iz,z') 



16^ 



d^k Tr,JGk(z) 



5^ 

dkj 

Gk(z)Gk(z')^(n) 



-2 ^ Gk(z) — 
dkf, dkf, 

Here Gk(z) = (//k - z) ' is the single-particle Greens function 
in reciprocal-space (plane-wave) representation. Application to 
graphene yields the same results as obtained more easily using 
formulae from the previous subsection. However, Eq. ^ turns 
out to be very useful beyond the clean limit. 

3. Single-layer graphene in the clean limit 

We apply the continuum-model description of a single sheet 
of graphene to evaluate the conductivity formula (|7]). The 
single-particle Hamiltonian in plane-wave representation is 
given by SH 

= hv {k_,o-_, + kyo-y) + T [{^kj - kf) cr, + Ik.kyO-y] , (12) 

where v is the Dirac-fermion velocity characterizing the dis- 
persion near the K point, and the term proportional to t is a 
trigonal-warping correction to the band structure. Straightfor- 
ward diagonalization of H^g yields eigenvalues e^^, where cr = 
+ distinguishes the electron and hole bands. As e.*''^' = -s't^^ , 

kcr k,-(T' 

the dependence on temperature and chemical potential is uni- 
versal, i.e., independent of the values of v and t. The remain- 
ing dimensionless prefactor is only a function of w, v, t and, 
by simple dimensional analysis, can therefore only depend on 
these through the combination ojt/ (Hv^). The expression for the 
conductivity reads then 
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Figure 1: Universal function capturing tlie effect of trigonal warping on tlie 
optical conductivity of a single graphene sheet. The dashed lines indicate the 
values of n/S and respectively. 



where is the Heaviside step function, and we introduced the 
abbreviation 



sinh^ 



cosh^ + cosh 77 
It is found that T^,, - T^, , - and 



(14) 



J 1 (l+2^^,)2(^,-^^2) 



X(5 (1 - 2 _ 2^K, [kI - 3/f2] + ^Vj . (15) 

We find numerically that T.„(^) = T^./^) = T(^) and show 
this universal function in Fig. [1] In the limit r — > 0, the 
well-known (l3 , 2^ 23] universal conductivity of single-layer 
graphene is found: T(0) - n/^. In the (for graphene unphys- 
ical) limit of large ^, a different saturation value is realised: 

T(oo) = 7T/4. 



4. Bilayer graphene in the clean limit 

To describe bilayer graphene, we use the 4x4 continuum- 
model Hamiltonian 1 16] in plane-wave representation given by 



Hhi 



eV/2 
hv(kx — iky) 
t± 




hv(kx + iky) 
eV/2 


Hv^ikx - iky) 









-ey/2 
hv{kx + iky) 



hv^ikji + iky) 
hv(kx — iky) 
-eV/2' , 
(16) 

Here and V parameterize the strongest inter-layer hopping 
and a potential difference applied between the layers, respec- 
tively. V3 measures the strength of an additional inter-layer 
hopping that gives rise to trigonal warping. Straightforward di- 
agonalization of //bi yields the set of energy eigenvalues s^^, 
with cr — 1 , 2, 3, 4. These eigenvalues can be grouped into two 
pairs that add up to zero. Assuming without loss of general- 



ity that e. 
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we have e, 



(bl) 
kl 



-e^^ and 



This means that there are two contributions of 
the type encountered in the single-layer case. Depending on the 
situation, the remaining four contributions can be simplified as 




Figure 2: Interband contribution to the conductivity of clean bilayer graphene 
{ti_ = lOk^T), with chemical potential at the symmetry point and a finite bias 
voltage V = Ik^T/e between the layers (thick solid curve). The zero-bias case 
is shown as the thin solid curve. The dashed line indicates the value 7t/4. 



well. Here we just give the analytical result obtained for the 
case with zero inter-layer bias and trigonal warping neglected 



(y = and V3 = 0). We find cr] 
where 



.(inter) 



- cr- 



.(inter) 

yy 



= 0-1+0-2+ 0-3, 



2k^'k^ 



fui) + 2ti_ 



hcj + 1± 
hto — 2tj_ 

hoj — ti_ 



0(w) 



11 . U^\\^^.,jn)U 

(To 8 \ tUi)l [ \ 



0(w - 2t^lK) 
hio 2ji - 1^^ 



(17) 



+g 



2kBT' 2ksT 
hoj 2/i + fj^ 



2kKT 2knT 



(18) 



0-3 tj_ r dK / f± ^ 

— = —6{oj-t^lh)\ — ^ kr-f ' 7—7; - 

0-0 h K [ \2kYiT kuT 



+8 



2kBT' ksT 



+ K 



(19) 



This result generalizes a previous expression [23] obtained for 
the zero-temperature limit. In Fig. |2] we show the effect of a 
finite inter-layer bias V. 



5. Conclusions 

We have studied the ac conductivity of single and bilayer 
graphene. Analytical results were obtained for finite temper- 
ature and with trigonal warping included. Our expressions 
should be useful to facilitate detailed comparison with exper- 
iment and enable extraction of electronic-structure parameters 
from conductivity measurements. 



Acknowledgment 

JZB is supported by a postdoctoral fellowship grant from the 
Massey University Research Fund. 



3 



References 



[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. 

Dubonos, I. V. Giigorieva, A. A. Firsov, Science 306 (2004) 666. 
[2] K. S. Novoselov, D. Jiang, F. Scliedin, T. J. Booth, V. V. Khotkevich, S. V. 

Morozov, A. K. Geim, PNAS 102 (2005) 10451. 
[3] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, A. K. 

Geim, Rev. Mod. Phys. 81 (2009) 109. 
[4] M. Katsnelson, K. Novoselov, Sol. St. Comm. 143 (2007) 3. 
[5] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 1. Katsnelson, 

I. V. Grigorieva, S. V. Dubonos, A. A. Firsov, Nature (London) 438 (2005) 

197. 

[6] Y. Zhang, Y. Tan, H. L. Stormer, P Kim, Nature (London) 438 (2005) 
201. 

[7] M. 1. Katsnelson, Eur. Phys. J. B 51 (2006) 157. 

[8] V. P Gusynin, S. G. Sharapov, J. R Carbotte, Phys. Rev. Lett. 96 (2006) 
256802. 

[9] V. R Gusynin, S. G. Sharapov, Phys. Rev. B 73 (2006) 245411. 
[10] V. P Gusynin, S. G. Sharapov, J. R Carbotte, Phys. Rev Lett. 98 (2007) 
157402. 

[11] V. P Gusynin, S. G. Sharapov, J. P Carbotte, Phys. Rev. B 75 (2007) 
165407. 

[12] K. Ziegler, Phys. Rev. Lett. 97 (2006) 266802. 
[13] K. Ziegler, Phys. Rev. B 75 (2007) 233407. 

[14] N. M. R. Peres, R Guinea, A. H. C. Neto, Phys. Rev. B 73 (2006) 125411. 
[15] J. Nilsson, A. H. Castro Neto, F. Guinea, N. M. R. Peres, Phys. Rev. Lett. 
97 (2006) 266801. 

[16] J. Nilsson, A. H. C. Neto, F Guinea, N. M. R. Peres, Phys. Rev. B 78 
(2008) 045405. 

[17] D. S. L. Abergel, V. I. Fal'ko, Phys. Rev. B 75 (2007) 155430. 
[18] J. Cserti, Phys. Rev B 75 (2007) 033405. 

[19] J. Cserti, A. Csordas, G. David, Phys. Rev. Lett. 99 (2007) 066802. 
[20] L. A. Falkovsky, A. A. Varlamov, Eur Phys. J. B 56 (2007) 281. 
[21] 1. Snyman, C. W. J. Beenakker, Phys. Rev. B 75 (2007) 045322. 
[22] S. Ryu, C. Mudry, A. Furusaki, A. W. W. Ludwig, Phys. Rev. B 75 (2007) 
205344. 

[23] E. J. Nicol, J. R Carbotte, Phys. Rev. B 77 (2008) 155409. 
[24] T. Stauber, N. M. R. Peres, A. K. Geim, Phys. Rev. B 78 (2008) 085432. 
[25] N. M. R. Peres, T. Stauber, Int. J. Mod. Phys. B 22 (2008) 2529. 
[26] N. M. R. Peres, T. Stauber, A. H. Castro Neto, Europhys. Lett. 84 (2008) 
38002. 

[27] L. A. Falkovsky, arXiv:0806.3663 

[28] H. Min, A. H. MacDonald, Phys. Rev. Lett. 103 (2009) 067402. 

[29] F. Wang, Y. Zhang, C. Tian, C. Girit, A. Zettl, M. Crommie, Y. R. Shen, 

Science 320 (2008) 206. 
[30] R. R. Nair, R Blake, A. N. Giigorenko, K. S. Novoselov, T. J. Booth, 

T. Stauber, N. M. R. Peres, A. K. Geim, Science 320 (2008) 308. 
[31] Z. Q. Li, E. A. Henriksen, Z. Jiang, Z. Hao, M. C. Martin, P Kim, H. L. 

Stormer, D. N. Basov, Nat. Phys. 4 (2008) 532. 
[32] K. R Mak, M. Y. Sfeir, Y. Wu, C. H. Lui, J. A. Misewich, T. R Heinz, 

Phys. Rev. Lett. 101 (2008) 196405. 
[33] L. M. Zhang, Z. Q. Li, D. N. Basov, M. M. Fogler, Z. Hao, M. C. Martin, 

Phys. Rev. B 78 (2008) 235408. 
[34] Z. Q. Li, E. A. Henriksen, Z. Jiang, Z. Hao, M. C. Martin, P Kim, H. L. 

Stormer, D. N. Basov, Phys. Rev. Lett. 102 (2009) 037403. 
[35] O. Madelung, Introduction to Solid-State Theory, Springer, Berlin, 1978. 
[36] J. C. Slonczewski, P. R. Weiss, Band structure of graphite, Phys. Rev. 109 

(1958) 272. 

[37] R. Saito, G. Dresselhaus, M. S. Dresselhaus, Physical Properties of Car- 
bon Nanotubes, Imperial College, London, 1998. 



4 



